Anomalous scaling in depinning transitions 
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It is demonstrated that the renormalization group (RG) flows of depinning transitions do not depend 
on whether the driving force or the system velocity is kept constant. This allows for a comparison 
between RG results and corresponding self-organized critical models. However, close to the critical 
point, scaling functions cross over to forms that can have singular behavior not seen in equilibrium 
thermal phase transitions. These can be different for the constant force and constant velocity driving 
modes, leading to different apparent critical exponents. This is illustrated by comparing extremal 
dynamics for interface depinning with RG results, deriving the change in apparent exponents. Thus 
care has to be exercised in such comparisons. 

PACS numbers: 68.35.Rh, 05.65.+b, 64.60.Ht, 05.45.-a 



Driven dynamics of disordered systems have been stud- 
ied extensively from two different approaches. As the ex- 
ternal driving force is increased, a system undergoes a 
depinning transition from a macroscopically static state 
(with transient motion) to a moving steady state. This 
transition has been successfully described as a dynami- 
cal second order phase transition |^ , and analyzed within 
the framework of the renormalization group (RG), orig- 
inally for charge-density waves and subsequently for 
interfaces and a variety of other systems |^,||. This 
approach obtains the scaling behavior near the transition 
where the system starts to move. On the other hand, the 
concept of self-organized criticality (SOC), in the original 
sandpile model Q and descendants thereof Q has been 
used to obtain the behavior of systems that are forced to 
stay at the depinning transition. 

One would expect there to be connections between 
these two viewpoints Q]; indeed, it is possible to 
map automaton models for CDWs to the original SOC 
model 1^,^, and 1-1-1 dimensional moving interface to a 
slightly different SOC model [^. The connections be- 
tween the two have been exploited to obtain the dynam- 
ical critical exponent z for CDWs by borrowing from ex- 
act results for SOC sandpiles jll]. However, it has been 
pointed out that one has to be careful whether the 
depinning transition is approached with a time indepen- 
dent external driving force, tuned to its critical value, or 
with a time independent (infinitesimal) velocity. In the 
latter case, one envisages a feedback loop that adjusts the 
external force in a time-dependent way, so as to ensure 
that the rate at which the system moves (averaged over 
its entire spatial extent) is strictly time independent. It 
has been suggested Q that the critical behavior at the 
depinning transition could be different for these two driv- 
ing modes. As an extreme example of constant current 
driving, there is the class of "extremal models" p2|-p^, 
where at every (discrete) time step, there is activity only 
at one lattice site in the system. 

For constant force driving, there are fluctuations in 
the (spatially averaged) velocity of the system, while for 



constant velocity driving there are fluctuations in the ex- 
ternal force. In either case, at a non-zero driving rate, the 
fluctuations are small for sufflciently large systems, and 
thus the two driving mechanisms should be equivalent. 
However, the large system limit is problematic at the 
critical point. For instance, with constant force driving, 
the mean square fluctuations in the velocity of a large sys- 
tem are (Sv)'^ ~ v'^{^/L)'^, where v is the mean velocity, 
^ is the correlation length at velocity v (defined through 
the velocity autocorrelation function or other methods), 
L is the linear extent of the system, and d is its dimen- 
sionality. One would expect that holding the velocity 
constant by adjusting the force will make a qualitative 
difference to the dynamics when the velocity fluctuations 
are a significant fraction of the mean velocity, i.e. f ~ L. 
Heuristically, for a fixed v and L — > co, imposing a con- 
stant velocity does not affect the dynamics, since the to- 
tal 'activity' at every time step is vL"^, so that avalanches 
etc. are free to unfold with their own internal timescales. 
On the other hand, for fixed L as ?; — > 0, imposing a 
constant velocity 'chokes off' the dynamics: avalanches 
proceed sequentially, with internal dynamics determined 
by the velocity constraint. 

In this paper, we show that (to the extent that con- 
tinuum descriptions are appropriate) both constant force 
and constant velocity driving actually have the same scal- 
ing under the renormalization group. However, scaling 
functions for various quantities, which can have different 
behavior in the limits ^ << L and ^ >> L, are often dif- 
ferent in the ^ >> L regime for the two driving modes. 
It is standard for scaling functions in equilibrium critical 
phenomena to behave differently when their arguments 
are small and large, but as will be discussed later, dynam- 
ical critical points have even more freedom in how scaling 
functions can behave. It is this freedom that allows iden- 
tical RG fiows to still yield different scaling functions for 
the two driving modes in the » L regime. 

We demonstrate this by focusing on a particular sys- 
tem, that of a interface with short range internal elastic 
interactions that is driven through a disordered medium 
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under the influence of an external force. Due to the dis- 
order, different parts of the interface experience different 
random pinning forces as they move forward. Using the 
RG relationship between the interface velocity and the 
correlation length [D w ~ i^^^' , we see that the crossover 
to ^ > L occurs at w ~ i^^^. Note that extremal dy- 
namics for a d -|- 1 dimensional interface corresponds to 
V ^ L~'^, and z — C,<dui all dimensions, so that ^ >> L 
for large L. 

It has been pointed out that in extremal 

models for the motion of pinned interfaces, the rough- 
ness of the interface can scale anomalously. Thus for an 
1-1-1 dimensional interface in a system of size L in steady 
state lH^], if h{x,t) is the interface position as a func- 
tion of the transverse coordinate x and the time t, the 
roughness w{t) = [{{h{x, t + to) — h{t + to) — h{x, to) + 
h{to)}^)]^/'^ (all averages are spatial averages) has the 
scaling form |l^,|l8| 



w{t)^t^^L~'/'J 



(1) 



(The subscripts on /3 and z denote constant velocity driv- 
ing.) It can be shown [ p^Jl6| that the scaling function ip 
goes to a constant for large values of its argument. The 
explicit L dependence that remains is anomalous, con- 
trary to the normal expectation of a well defined L — > oo 
limit. As a consequence of this, if the long time (steady 
state) roughness scales as wit — > oo) ~ L^, then the 
conventional relation ^ = z/3 is replaced with 
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With the additional resuh ISUlflJl 



1 + C 



(2) 



(3) 



one can obtain [3y in terms of Q. Both Eq.(g) and Eq.(||) 
differ from renormalization group (RG) results [| for an 
interface driven with a constant applied force. In particu- 
lar, the RG analysis yields no simple relation between the 
dynamic exponent z and the roughness exponent Cj un- 
like Eq.(|^). Nor is there any reason to expect the anoma- 
lous scaling form of Eq.(|l|), and therefore C, = z(3 should 
hold instead of Eq. (||) . 

We first demonstrate that, despite appearances, the 
RG flows are unaffected by going from a constant force 
driving mode to a constant velocity driving mode. The 
standard equation for zero-temperature driven dynamics 
of an interface is 120(1 



dth{x, t) = V^h{x, t) + Y{h{x, t);x)+ F. 



(4) 



Here \/'^h{x, t) comes from the (short-ranged) elastic en- 
ergy of the interface, Y{h; x) is a pinning force that comes 
from a random (impurity) potential, and F is the external 



driving force. (Of course, it is not necessary that all lat- 
tice growth models — e.g. the Sneppen model [T^ ] — can 
be mapped to this or any continuum equation.) We have 
neglected KPZ like terms pl[. Using the Martin Sig- 



gia Rose method |23 and introducing an auxiliary field 
h(x,t), one can construct a generating functional Z 
which can be written as 



[dh] [dh] exp 



d'^xdt 



ih{x, t){dth - y^h ~F- Y{h; x)} 



(5) 



where integrating out the auxiliary field yields a prod- 
uct of i5-functions that imposes Eq.(^). If instead the 
interface is driven at constant velocity, Eq. (^) is replaced 
with 



Z' = J [dF] [dn] [dh] [dh] exp J d'^x dt 



ih(x, t){dth - V^/i - F{t) - Y{h; x)} 
+ifi{t){v — dth}. 



(6) 



This extension is actually not difficult to understand: at 
any time t, by first integrating over /i(i) and h(x,t), we 
obtain (5- function constraints that impose (dth) — v (the 
average here is a spatial average) and Eq.(|^) with a (as 
yet unknown) driving force F(t). Now performing the 
integral over F, the integral together with the first con- 
straint sets F{t) to be whatever it has to be for {dth to 
be equal to v. One is left with the second constraint, i.e. 
Eq.(^, with F{t) adjusted to ensure constant velocity. 
Even though F{t) is now a dynamical variable instead of 
a parameter, since it is not a field, i.e. it has no spatial 
dependence (nor has fi(t)), when short distance fluctua- 
tions are eliminated under renormalization there are no 
extra loop corrections in Z' compared to Z. 

As mentioned earlier, even though the RG flxed point 
is the same for constant force and constant velocity driv- 
ing, the behavior of scaling functions are different for the 
two driving modes in the ^ >> L regime. To illustrate 
this, we first consider the scaling of the duration of an 
avalanche as a function of its linear size. This has the 
general scaling form 



t[lLX) = l'T{l/L,^/L). 



(7) 



Here / is the linear extent of the avalanche (in the d trans- 
verse directions; in the direction the interface moves, the 
extent is ~ K), L is the linear size of the system, ^ is 
the correlation length, and z is the non-trivial dynamical 
exponent from the RG. For a fixed I and ^, as L — + oo, 
the avalanche duration should be independent of L since 
the system is uncorrelated on length scales much bigger 
than ^ It should also be independent of ^, since 

requiring a total growth rate of vL'^ for the system does 
not affect any individual avalanche from proceeding with 
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its own intrinsic timescale. Thus r(0, 0) is some non-zero 
constant. 

In the other regime of ^ >> L, i.e. if the u ^ limit is 
taken before L — > oo, the avalanches are non-overlapping 
in time. In the constant velocity driving mode, where the 
velocity constraint is imposed, any single avalanche pro- 
ceeds at a fixed rate. Therefore for the constant velocity 
driving mode, til, L, ^) must be inversely proportional to 
the rate at which the avalanche is allowed to proceed, 
which is vL"^. Requiring that t should have such a de- 
pendence on V and L, using the result [|| w ~ it is 
straightforward to verify that 



with 



lim t{l,L,i) 

^/L— >oo 



l^+'^l{vL'^). 



(8) 



This yields an apparent dynamical exponent of z^, = C,-\-d, 
which is the d -\- 1 dimensional generalization of Eq. . 
Note that there is no change in the apparent dynamical 
exponent in the ^ » L regime for constant force driving, 
where an avalanche is allowed to proceed at its intrinsic 
rate (through parallel updating of lattice sites). This is 
why the RG estimate for z for two dimensional CDWs 
agrees with the exact result for Abelian sandpiles . 

The analysis of the interface roughness for a system of 
linear size L, w{t, L), proceeds in a similar manner. At a 
velocity v, the interface roughness has the scaling form 



w{t,L,S,)^L^W{t/L'--i/L). 



(9) 



For fixed t and ^, the roughness must have a well defined 
L ^ oo limit: 



(10) 



For large t, the roughness saturates to the steady state 
form w ^ ^'^ . For ^ ^ cxd (the large L limit has been 
taken first), or equivalently for small the roughness is 
^ independent, i.e. w ~ t^^^ . 

In the other regime of f >> L, the apparent exponents 
are once again different for constant velocity driving. In 
this case, the dynamics are controlled by the imposed 
velocity, i.e. the v and t dependence of Eq.(H) occurs 
only in the combination r = vL'^t. Using v ~ this 
implies 



lim ' 



•j{t,L,i)^L^W2{T/L^^ 



(11) 



(With extremal dynamics, where one site moves forward 
at every time step, r — t.) In the large time limit, the 
roughness approaches the steady state form w ^ L'' , i.e. 
W2{oo) is a constant. It is not obvious how to extract 
the behavior of the function W2 when its argument is 
small, but physical arguments supported by numerical 
results [|l3|,|l6| show that w{t, L) must have a residual 
]^-d/2 dependence. Therefore 



Pv = 



C + d/2 

C + d 



(13) 



which is the same as Eq.(lO) of Ref. (generalized to 
d dimensions). 

Eqs.(||) and ( |l^ ) give the w ^ scaling behavior of 
avalanche durations and interface roughness respectively 
for constant velocity driving. We note once again that 
these are for ^ >> L, i.e. w — > before L — » 00. As men- 
tioned before, even in critical phenomena for equilibrium 
phase transitions, as one approaches the transition for a 
fixed system size, one sees a change in the scaling form 
of dynamical variables. However, it is generally possi- 
ble to obtain the behavior in this regime by requiring 
that there should be no dependence on (say) the reduced 
temperature t for the behavior of a finite size system. 
This requirement comes from the fact that there are no 
thermodynamic singularities for a finite sized system: it 
is possible to go smoothly from one side of the phase 
transition to another. No such requirement exists for dy- 
namical phase transitions, and one must be careful about 
possible u-dependence even in the u — > regime [||. 

There are other mechanisms as well that can make the 
connection between RG exponents and apparent scaling 
difficult. For instance, for CDWs below the depinning 
threshold, the periodicity of the dynamical variable (the 
CDW phase) prevents it from advancing by more than 27r 
anywhere in a single avalanche, but the same periodicity 
makes a region that has just avalanched susceptible to an 
imminent 'retriggering' of a fresh avalanche The low 
frequency dynamics thus sees a non-trivial analog of the 
roughness exponent, and the distribution of avalanche 
sizes is singular as u — > 0. Another example is for interface 
roughness itself, where the steady state roughness over a 
subsystem of size a; in a system with size L scales as 
w'^{x;L) ~ x^L^^"^ when C > 1 lUl- This is because 
the Fourier transform of the steady state roughness must 
be well behaved as L — + cx), and for C > 1 w'^{x;L) is 
dominated by long wavelength modes with q 1/L. 

In this paper we have shown that the renormaliza- 
tion group fiows for depinning transitions are the same 
whether the system is driven with a constant force or a 
constant velocity. However, in the critical regime, the 
apparent scaling behavior of physical quantities can be 
different for the two driving modes. This is because scal- 
ing functions can depend on the system velocity in a sin- 
gular manner as the transition is approached. This is 
unlike what is seen in equilibrium critical phenomena, 
and seems to be more common in the constant velocity 
driving mode. Despite this caveat, it is possible to obtain 
results for one driving mode from the other. 

I thank Anne Tanguy and Maya Paczuski for useful 
discussions. 
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